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Abstract 

SPYTHIA is an event level Monte Carlo program which simulates particle production 
and decay at lepton and hadron colliders in the Minimal Supersymmetric Standard 
Model (MSSM). It is an extension of PYTHIA 5.7, with all of its previous capabilities. 
This paper is meant to supplement the PYTHIA/ JETSET user manual, providing a de- 
scription of the new particle spectrum, hard scattering processes, and decay modes. 
Several examples of using the program are provided. 
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1 Introduction 



Recently, experiments at the Tevatron collider at Fermilab and the Large Electron Positron 
(LEP) collider at CERN have greatly extended our understanding of the Standard Model 
(SM). Several LEP energy upgrades have been staged on a migration towards LEP190, which 
will be replaced by the Large Hadron Collider (LHC). Prior to the operation of the LHC, the 
Tevatron collider will begin the Main Injector era, possibly with an extended running time, 
perhaps even at a higher luminosity than originally planned. Despite its phenomenological 
successes, the SM suffers from various theoretical problems which make it seem unlikely to 
be a complete theory. The rich physics program outline above is expected to definitively 
probe the energy regime responsible for Electroweak Symmetry Breaking (EWSB). One 
promising extension of the SM, which successfully addresses several of its deficiencies, is the 
the Minimal Supersymmetric Standard Model (MSSM) Q. Monte Carlo simulations are 
useful tools for studying the phenomenological implications of a new theory, particularly 
for determining search strategies and optimizing detector design. SPYTHIA is an event 
level Monte Carlo program which simulates particle production and decay at lepton and 
hadron colliders in the MSSM. It is an extension of PYTHIA 5.7 ||, with all of its previous 
capabilities, and additional particles, hard scattering processes, and decays. Furthermore, 
the simulation of the MSSM Higgs sector already present in PYTHIA is extended to include 
the decay of Higgs bosons into MSSM particles and the decay of MSSM particles into Higgs 
boson. 

The are already a few similar programs available. Most of the processes included in 
SPYTHIA are already present in the ISAJET/ISASUSY || event generator. In addition, a 
number of processes and decays relevant for lepton colliders are included in the generator 
SUSYGEN 0, which is interfaced to JETSET || and includes initial state photon radiation. 
The development of SPYTHIA arose from the desire to study the phenomenology of the 
MSSM at lepton-lepton, lepton-hadron, and hadron-hadron colliders using the initial and 
final state radiation and fragmentation models of PYTHIA/ JETSET. Additionally, one can 
perform cross checks and estimate model uncertainties with other generators. 

Like the SM, the MSSM contains a number of parameters with the dimension of mass 
which are not fixed by the theory. Supergravity (SUGRA) inspired models reduce the 
number of free parameters by imposing universality and exploiting the apparent unifica- 
tion of gauge couplings. Five parameters fixed at the gauge coupling unification scale, 
tan/3, Mo, mi/2i an d sign(n), are then related to the mass parameters at the scale of 
EWSB by renormalization group equations ||. ISASUSY and SUSYGEN numerically solve 
these equations to determine the mass parameters. Alternatively, they can input a general 
set of parameters. SPYTHIA operates in the second manner, with a slightly more general 
set of input parameters. There are three reasons for this: (1) programs already exist to 
calculate the SUGRA inspired mass parameters, (2) approximate analytic formulae |jj also 
exist which reproduce the output of ISASUSY within ~ 10%, and (3) we desire to study a 
much richer phenomenology than that possible in SUGRA inspired models. The SPYTHIA 
input parameters are described in detail later. 

In the following, it is assumed that the reader has a working knowledge of PYTHIA/ JETSET^. 
The modifications made in SPYTHIA are the sole product and responsibility of this author 
and not H.U. Bengtsson or T. Sjostrand. However, these modifications will be included 

2 A postscript version of the user manual f or PYTHIA 5. 7/ JETSET 7.4 is available at the URL address 
http: / /thep. lu.se /tf2 / staff/torbjorn / manual.ps . 
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in the future release of PYTHIA 6.1. Sec. 2 explains the physics assumptions behind the 
implementation of supersymmetry in PYTHIA and catalogs the new particles, processes and 
input parameters. The parameters and routines used in the MSSM simulation are described 
in Sec. 3. Availability and setup of the SPYTHIA package are detailed in Sec. 4, along with 
some examples of its use. Conclusions are presented in Sec. 5. 



2 Simulation of Supersymmetry 
2.1 Particle Spectrum 

SPYTHIA assumes the particle content of the MSSM. Each SM fermion has a scalar partner 
with the same quantum numbers. For counting purposes, the fermion states tpL and ipR 
are separate particles with scalar partners 0^ and Each SM gauge boson has a fermion 
partner with the same quantum numbers. The Higgs sector is extended to two complex 
scalar doublets, leading to 3 additional physical Higgs bosons H, A and H 1 * 1 to complement 
the SM Higgs h. Fermion higgsino partners are added for the two scalar doublets. Addi- 
tionally, a light gravitino G is included to allow studies of gauge mediated dynamical SUSY 
breaking ||]. The particle partners and SPYTHIA KF codes are listed in Table 1. Note that 
antiparticles of scalar particles are denoted by * in the text. 

Because of the large Yukawa couplings and the running masses for colored particles, 
the interaction and mass eigenstates for the third generation sfermions can be significantly 
mixed. We denote the top and bottom squark mass eigenstates ti,t2,bi, and b2 to distinguish 
them from the nearly degenerate eigenstates for the lighter squarks. In addition, the tau 
slepton mass eigenstates are f\ and r-i- The subscript 1 or 2 refers to the lightest and heaviest 
state respectively. In SUGRA inspired models, in the absence of mixing so that interaction 
eigenstates are the same as mass eigenstates, the right eigenstate is lighter than the left. In 
this case, for example, bi = b# and hi = b^. For completeness, we include the sterile Dr 
particles. The mixing of the partners to the electroweak gauge bosons (gauginos) and the 
two Higgs doublets (higgsinos) lead to the mass eigenstates for neutralino Nj and chargino 
Cj particles. A complete listing of the full particle spectrum, including decay channels, 
branching ratios, and the total width is available (as usual) through the subroutine LULIST. 



2.2 Particle Decays 

In PYTHIA, resonances are decayed so that color flow information is correctly passed from 
unstable particles to their decay products and no unstable colored particles are passed on 
to JETSET. In SPYTHIA, all Supersymmetric partners are treated as resonances. SPYTHIA 
assumes R-symmetry, which has two major phenomenological consequences: (1) Supersym- 
metric particles are produced in pairs, and (2) there is one and only one stable Supersym- 
metric particle. By default, the lightest superpartner (LSP) is Ni, unless Ni is allowed to 
decay to a gravitino. In that case, Ni is the next to lightest superpartner (NLSP) and G is 
the LSP. 



The decays of superpartners are calculated using the formulae of Refs. |10,|il, 12, 14]. All 
decays are spin averaged. For simplicity, the kinematics of three body decays of neutralinos, 
charginos, and the gluino are sampled using only the phase space weight. Decays involving 
b and t use the formulae of O], so they are valid for large values of tan/3. The one loop 
decays Nj — > Nj7 and t — » cNi are also included. 
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Table 1: SPYTHIA particle KF codes. 



KF 


Name 


Printed 


KF 


Name 


Printed 


A 1 

41 


1 


~d_L 


59 




~numu_L 


42 


1 

d« 


~d_R 


60 


V/J.R 


~numu_R 


A O 
4:6 


UL 


~u_L 


61 


T\ 


~tau_l 


A A 

44 


Ui? 


~u_R 


OZ 


T2 


~tau_2 


A FC 


§1 


_ T 

S_L 


DO 


V T L 


nutau_L 


4±0 


S2 


S_n, 


04fc 


v tR 


ITUX a.U_Tl 


47 


L L 


~r T 






g 


A Q 

4o 




C_H 


00 


AT 

JNi 


Ni 


A n 

49 


bi 


~b_l 


67 


AT 

IN 2 


~N 2 


50 




~b_2 


68 


AT 

JN 3 


~N 3 


k 1 
ol 


J- 

tl 


~4_ -\ 
t_l 


fin 
09 


AT 

IN 4 


JM4 


52 


t 2 


~t_2 


70 


cf 


"Ci 


53 


ei 


~e_L 


71 




~c 2 


54 




~e_R 








55 




~nue_L 








56 




~nue_R 








57 


Al 


~mu_L 








58 


VR 


~mu_R 


76 


G 


"grvtno 



2.3 Processes 

Tables 2-5 are meant to update Tables 11-14 in the Pythia 5.7 manual. The number of 
available processes in SPYTHIA is extended to 400, with processes 201-280 reserved for the 
MSSM. In processes 210 and 213, t refers to both e and jl. For ease of readability, we 
have removed the subscript L on u. tjt*,fjT* and fii 1 * production correctly account for 
sfermion mixing. Several processes are conspicuously absent from the table. For example, 
processes 255 and 257 would simulate the associated production of right handed squarks 
with charginos. Since the right handed squark only couples to the higgsino component of 
the chargino, the interaction strength is proportional to the quark mass, so these processes 
can be ignored. 

Because there are so many processes involved, there are shortcuts to allow the simulation 
of various classes of signals. These classes are accessible through the parameter MSEL, and 
are listed in Table 6. 
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Table 2: Subprocess codes, part 1. First column is '+' for processes implemented and blank 
for those that are only foreseen. Second is the subprocess number ISUB, and third the 
description of the process. The final column gives references from which the cross sections 
have been obtained. 



In 



No. 



Subprocess 



Reference 



+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 



201 
202 
203 
204 
205 
206 
207 
208 
209 
210 
211 
212 
213 
214 
216 
217 
218 
219 
220 
221 
222 
223 
224 
225 
226 
227 
228 
229 
230 
231 
232 
233 
234 
235 
236 



UU 
UU 

fif; 
uu 
uu 

UU 

uu 
uu 
uu 

f-f • 
f-f • 

tfi 
uu 
uu 
uu 
uu 
uu 
uu 
uu 
uu 
uu 
uu 
uu 
uu 
uu 
uu 
uu 

f-f • 
f-f • 
f-f • 
f-f • 
f-f • 
f-f • 
f-f • 
f-f • 



e L e* R + e* L e R 

T2T 2 

■ hv} + ? L Vl 

■ T 2 l>r* + T%V T 

~ ~ * 

v T v T 

NiNi 

N 2 N 2 

N 3 N 3 

N4N4 

NiN 2 

N1N3 

N1N4 

N 2 N 3 

N 2 N 4 

N3N4 

cfcf 

p±p=F 

cfc^ 

-NiC± 
■N 2 Cf 
N 3 Cf 
-N 4 C± 
■NxCf 

- N 4 Cf 



B 

13, 

Hi 



[13 



m 




El 

15 



m 
p 



q 

Pi 





El 




















El 




E 



E 
E 



E 
E 

E 

11 
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Table 3: Subprocess codes, part 2. First column is '+' for processes implemented and blank 
for those that are only foreseen. Second is the subprocess number ISUB, and third the 
description of the process. The final column gives references from which the cross sections 
have been obtained. 



In 



No. 



Subprocess 



Reference 



+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 



237 
238 
239 
240 
241 
242 
243 
244 
246 
247 
248 
249 
250 
251 
252 
253 
254 
256 
258 
259 
261 
262 
263 
264 
265 
271 
272 
273 
274 
275 
276 
277 
278 
279 
280 



kU -» gNi 
:N 2 
N 3 

N 4 
± 



% 
uu 



ff 
ff • 

fif,: 



fig 
fig 
fig 
fig 
fig 
fig 
fig 
fig 
fig 



gCf 

gcf 



QiLNi 

qifiNi 

QiLN 2 

QiLN 3 
qii?N 3 



f*g-> 




Sg-> 


qiLg 


Hg^ 


CLiRg 


uu- 


'hi* 


uu- 


> t 2 t 2 


uu- 


► tlt 2 + tjt 2 


gg -» 


tltj 


gg -> 


t 2 t 2 


tfi - 


► qiLQjL 


Utj - 


► QiRCijR 


ff • — 


► QiLQ_jR + q«i?qjL 


ff — 


> qiL<ljL 


ff • — 


* qiiiq^i? 


% ~ 


* <iiL<ljR + Q.iRQLjL 


uu- 




uu- 


" QjRQjR 


gg 


q«Lq*L 


gg -* 


qi«q*-R 



© 
© 



© 
© 



HI 

15 
15 

HI 



© 
© 



© 
© 



© 
© 



© 

PI 



[15] 
PI 



© 
15 
15 



© 
PI 
15 



15 
15 
PI 
= 

mi 
mi 
mi 
mi 
mi 

15 
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Table 4: Classes of processes accessible through the parameter MSEL and the individual 
processes codes I SUB. 



MSEL 


Description 


39 


All MSSM processes except Higgs production 






40 


Squark and gluino production, ISUB = 243, 244, 258, 259, 


271- 


-280 


41 


Stop pair production, ISUB = 261-265 






42 


Slepton pair production, ISUB = 201-214 






43 


Squark or gluino with chargino or neutralino, ISUB = 237- 


-242, 


246-256 


44 


Chargino-neutralino pair production, ISUB = 216-236 







3 The Parameters and Routines of the MSSM Simulation 

All of the Supersymmetric extensions of the code are included in the standard PYTHIA 
library. By default, Supersymmetry is not simulated. However, by setting various parame- 
ters, a rich MSSM phenomenology is available. 

The parameters available to the user are stored in the FORTRAN common block 
C0MM0N/PYMSSM/IMSS(O:99) ,RMSS(0:99). In general, options are set by the IMSS array, 
while real valued parameters are set by RMSS. The entries IMSS(O) and RMSS(O) are not 
used, but are available for compatibility with the C programming language. The arrays are 
described below. The default values are denoted by (D). 



3.1 The MSSM Parameters 



COMMON/PYMSSM/IMSS (0 : 99) , RMSS (0 : 99) 



Purpose: to give access to parameters that allow the simulation of the MSSM. 
IMSS(l) : (D=0) level of MSSM simulation. 
= : No MSSM simulation. 

= 1 : A general MSSM simulation. The parameters of the model are set by the 
array RMSS. 

= 2 : An approximate SUGRA simulation using the analytic formulate of (7j 

to reduce the number of free parameters. In this case, only five input 
parameters are used. RMSS(l) is the common gaugino mass mi/2, RMSS (8) 
is the common scalar mass too, RMSS (4) fixes the sign of the higgsino 
mass (J,, RMSS (16) is the common trilinear coupling A, and RMSS (5) is 
tan/3 = vijv\. 

IMSS (2) : (D=0) treatment of C/(l), SU(2), and SU(3) gaugino mass parameters. 

= : The gaugino parameters M l7 M 2 and M 3 are set by RMSS(l) , RMSS (2) , 

and RMSS (3), i.e. there is no forced relation between them. 

= 1 : The gaugino parameters are fixed by the relation Mi/«i3/5 = M<ija<i = 
M 3 /a 3 = X and the parameter RMSS(l). If IMSS (1) =2, then RMSS(l) 
is treated as the common gaugino mass m X j 2 and RMSS (20) is the GUT 
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= 2 

IMSSC3) : 
= 



= 1 

IMSS(4) : 
= 

= 1 



= 2 
IMSS(7) : 

= 
= 1 



IMSS(8) : 
= 
= 1 

IMSS(9) : 

= 

= 1 
IMSS(IO) 
= 

= 1 

IMSS(ll) 
= 
= 1 

RMSS(l) : 



scale coupling constant acuTi so that X = mi/ 2 / acuT- 
Mi is set by RMSS(l), M 2 by RMSS(2) and M 3 = M 2 a 3 /a 2 . In such a 
scenario, the U(l) gaugino mass behaves anomalously. 
(D=0) treatment of the gluino mass parameter. 

The gluino mass parameter M3 is used to calculate the gluino pole mass 



(D= 



with the formulae of [16]. The effects of squark loops can significantly 
shift the mass. 

M3 is the gluino pole mass. The effects of squark loops are assumed to 
have been included in this value. 
1) treatment of the Higgs sector. 
The Higgs sector is determined by the approximate formulae of || and 
the pseudoscalar mass M4 set by RMSS(19). 

The Higgs sector is determined by the exact formulae of || and the pseu- 
doscalar mass Ma set by RMSS(19). The pole mass for Ma is not the 
same as the input parameter . 

The Higgs sector is fixed by the mixing angle a set by RMSS(18) and the 

mass values PMAS(I,1), where 1=25,35,36, and 37. 
(D=0) treatment of the scalar masses in an extension of SUGRA models. The 
presence of additional £7(1) symmetries at high energy scales can modify the 
boundary conditions for the scalar masses at the unification scale. 

No additional U-terms are included. In SUGRA models, all scalars have 

the mass mo at the unification scale. 

RMSS(23 — 25) are the values of Dx,Dy and Ds at the unification scale 

in the model of [p^ j. The boundary conditions for the scalar masses 

are shifted based on their quantum numbers under the additional £7(1) 

symmetries. 
(D=l) treatment of the f mass eigenstates. 

The f mass eigenstates are calculated using the parameters RMSS ( 13 , 14 , 17) . 

The r mass eigenstates are identical to the interaction eigenstates, so they 

are treated identically to e and fi . 
(D=0) treatment of the right handed squark mass eigenstates for the first two 
generations. 

The qR masses are fixed by RMSS (9) . d/j and ur are identical except for 
Electroweak D-term contributions. 

The masses of d# and ur are fixed by RMSS (9) and RMSS (22) respectively. 
(D=0) allowed decays for N2. 

The second lightest neutralino N2 decays with a branching ratio calculated 
from the MSSM parameters. 

N2 is forced to decay only to N17, regardless of the actual branching ratio. 
This can be used for detailed studies of this particular final state. 
(D=0) choice of the lightest superpartner (LSP). 
Ni is the LSP. 

Ni is the next to lightest superparter (NLSP) and the gravitino is the 
LSP. The gravitino decay length is calculated from the gravitino mass set 
by RMSS (21) and the Ni mass and mixings. 

If IMSS(1)=1 Mi, the U(l) gaugino mass. If IMSS(1)=2, then the common 

gaugino mass m^. 
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RMSS(2) 
RMSS(3) 
RMSS(4) 
RMSS(5) 
RMSS(6) 
RMSS(7) 
RMSS(8) 
RMSSC9) 
RMSS(IO) 
RMSS(ll) 
RMSS (12) 
RMSS (13) 
RMSS (14) 
RMSS (15) 
RMSS (16) 
RMSS (17) 
RMSS (18) 

RMSS (19) 
RMSS (20) 
RMSS (21) 
RMSS (22) 
RMSS (23) 
RMSS (24) 
RMSS (25) 



M2, the SU(2) gaugino mass. 

M3, the SU(3) (gluino) mass parameter. 

fi, the higgsino mass parameter. If IMSS(1)=2, only the sign of fi is used, 
tan j3, the ratio of Higgs expectation values. 

Left slepton mass Mg . The sneutrino mass is fixed by a sum rule. 
Right slepton mass . 

Left squark mass Mq L . If IMSS(1)=2, the common scalar mass mo- 
Right squark mass M^ R . Mq R when IMSS(9)=1. 
Left squark mass for the third generation Mq L . 
Right sbottom mass Mr . 
Right stop mass M~ tR 



Left stau mass M 



If negative, then it is assumed that M~ 2 ^ < 0. 



1~L • 



Right stau mass Mf R . 
Bottom trilinear coupling A^. 

Top trilinear coupling A t . If IMSS(1)=2, the common trilinear coupling A. 
Tau trilinear coupling A T . 

Higgs mixing angle a. This is only used when all of the Higgs parameters are 
set by the user, i.e IMSS(4)=2. 
: Pseudoscalar Higgs mass parameter M\. 
: (D=.041) GUT scale coupling constant cigut- 
: (D=1.0) The gravitino mass (eV). 
: ur mass when IMSS(9)=1. 

: Dx contribution to scalar masses when IMSS(7)=1 (GeV 2 ). 
: Dy contribution to scalar masses when IMSS(7)=1 (GeV 2 ). 
: Dg contribution to scalar masses when IMSS(7)=1 (GeV 2 ). 



C0MM0N/PYSSMT/ZMIX(4,4) ,UMIX(2,2) ,VMIX(2,2) ,SMZ(4) ,SMW(2) ,SFMIX(16,4) 



Purpose: to provide information on the neutralino, chargino, and sfermion mixing pa- 
rameters. The variables should not be changed by the user. 

ZMIX(4,4) : the neutralino mixing matrix in the Bino-neutral Wino-Up higgsino-Down 
higgsino basis. 

UMIX(2,2) : the chargino mixing matrix in the charged Wino-charged higgsino basis. 
VMIX(2,2) : the charged conjugate chargino mixing matrix in the wino-charged higgsino 
basis. 

SMZ(4) : the signed masses of the neutralinos. 
SMW(2) : the signed masses of the charginos. 

SFMIX(16,4) : the sfermion mixing matrices T in the L-R basis, identified by the corre- 
sponding fermion, i.e. SFMIX(6,I) is the stop mixing matrix. The four entries 
for each sfermion are Tn,Ti2,T2i, and T22- 



3.2 The MSSM Physics Routines 

The following subroutines and functions need not be accessed by the user, but are described 
for completeness. 

SUBROUTINE PYAPPS : uses approximate analytic formulae to determine the full set of 
MSSM parameters from SUGRA inputs. 
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SUBROUTINE PYGLUI : calculates gluino decay modes. 

SUBROUTINE PYGQQB : calculates three body decays of gluinos into neutralinos or charginos 
and third generation fermions. These routines are valid for large values of tan (3. 
SUBROUTINE PYCJDC : calculates the chargino decay modes. 

SUBROUTINE PYHEXT : calculates the non-Standard Model decay modes of the Higgs bosons. 
SUBROUTINE PYHGGM : determines the Higgs boson mass spectrum using several inputs. 
SUBROUTINE PYINOM : finds the mass eigenstates and mixing matrices for the charginos 

and neutralinos. 
SUBROUTINE PYMSIN : initializes the MSSM simulation. 
SUBROUTINE PYNJDC : calculates neutralino decay modes. 

SUBROUTINE PYPOLE : computes the Higgs boson masses using a renormalization group 

improved leading-log approximation and two loop leading-log corrections. 
SUBROUTINE PYRNMT : determines the running mass of the top quark. 
SUBROUTINE PYSFDC : calculates sfermion decay modes. 

SUBROUTINE PYSUBH : computes the Higgs boson masses using only renormalization group 
improved formulae. 

SUBROUTINE PYTBDY : samples the phase space for three body decays of neutralinos, charginos, 
and the gluino. 

SUBROUTINE PYTHRG : computes the masses and mixing matrices of the third generation 
sfermions. 

4 Setup and Use of the SPYTHIA Generator 

The SPTYHIA code is available as a uuencoded file on the World Wide Web0. It operates 
in the same fashion as the standard PYTHIA distribution. See the PYTHIA manual for a 
detailed example of its usage. Below, we present several examples of how to simulate 
MSSM phenomenology with SPYTHIA. 

4.1 Example 1: Light Stop 

The first example is an MSSM model with a light neutralino Ni and a light stop ti, so that 
t — > tiNi can occur. The input parameters are IMSS(1)=1, RMSS(1)=70., RMSS(2)=70., 
RMSS (3) =225 . , RMSS (4) =-40 . , RMSS (5) =1 . 5, RMSS (6) =100 . , RMSS (7) =125 . , RMSS (8) =250 . , 
RMSS(9)=250., RMSS (10) =1500., RMSS (11) =1500. , RMSS (12) =-128. , RMSS (13) =100 . , 
RMSS (14) =125., RMSS (15) =800., RMSS (16) =800 . , RMSS ( 17) =0., and RMSS (19) =400.0. 
The top mass is fixed at 175 GeV, PMAS(6 , 1)=175 . 0. The resulting model has = 55 
GeV and = 38 GeV. IMSS(1)=1 turns on the MSSM simulation. By default, there are 
no intrinsic relations between the gaugino masses, so M\ = 70 GeV, Mi = 70 GeV, and 
M3 = 225 GeV. The pole mass of the gluino is slightly higher than the parameter M3, and 
the decay g — > t*t + tit occurs almost 100% of the time. 

4.2 Example 2: Approximate SUGRA 

The second example is an approximate SUGRA model. The input parameters are IMSS ( 1 ) =2, 
RMSS(1)=200., RMSS (4) =1., RMSS (5) =10., RMSS(8)=800 . , and RMSS(16)=0.0. The re- 
sulting model has M- Al = 901 GeV, M &R = 890 GeV, M~ ti = 538 GeV, Mg L = 814 GeV, 

3 http: //www. hep. anl.gov/theory/mrenna/spythia.hmtl 
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Mg = 560 GeV, = 80 GeV, M e ± = 151 GeV, M h = 110 GeV, and M A = 883 GeV. 
It corresponds to the choice M =800 GeV, M 1/2 =200 GeV, tan/? = 10, A = 0, and 
sign(fi) > 0. The output is similar to an ISASUSY run, but there is not exact agreement. 

4.3 Example 3: ISASUSY Model 

The final example demonstrates how to convert the output of an ISASUSY^ run using the 
same SUGRA inputs into the SPYTHIA format. This requires a general model with the 
parameters IMSS(1)=1, IMSS(3)=1, RMSS(1)=83.81, RMSS(2)=168.90, RMSS(3)=581 .83, 
RMSS(4)=283.37, RMSS(5)=10., RMSS (6) =813 . 63, RMSS (7) =804 . 87, RMSS (8) =917 . 73, 
RMSS(9)=909.89,RMSS(10)=772.87,RMSS(11)=901.52,RMSS(12)=588.33,RMSS(13)=813.63, 
RMSS (14) =804. 87, RMSS (15) =610 . 54, RMSS(16)=422.35, RMSS (17) =600 . , and 
RMSS(19)=858.412. 



5 Conclusions 

SPYTHIA simulates the physics of the MSSM using the PYTHIA/ JETSET platform. References 
to the underlying physics that is simulated have been provided, as well as an overview of 
the user interface to the program. With the appropriate choice of parameters, a wide range 
of MSSM models can be studied, including those with a light gravitino and extra D-term 
contributions to scalar masses. 



information on the ISAJET and ISASUSY p rograms can be found on the World Wide Web at 
http://wwwcnl.cern.ch/asd/cernlib/mc/isajet.html[ 
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